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UNDERSTANDING GALAXY OUTFLOWS AS THE PRODUCT OF UNSTABLE TURBULENT SUPPORT 
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ABSTRACT 

The interstellar medium is a multiphase gas in which turbulent support is as important as thermal 
pressure. Sustaining this configuration requires both continuous turbulent stirring and continuous 
radiative cooling to match the decay of turbulent energy. While this equilibrium can persist for small 
turbulent velocities, if the one-dimensional velocity dispersion is larger than « 35 km/s, the gas moves 
into an unstable regime that leads to rapid heating. I study the implications of this turbulent runaway, 
showing that it causes a hot gas outflow to form in all galaxies with a gas surface density above 50 
Mq pc~^, corresponding to a star formation rate per unit area of 0.1 M0 yr~^ kpc~^. For galaxies 



with fo 



> 



200 km/s, the sonic point of this hot outflow should lie interior to the region containing 



cold gas and stars, while for galaxies with smaller escape velocities, the sonic point should lie outside 
this region. This leads to efficient cold cloud acceleration in higher mass galaxies, while in lower mass 
galaxies, clouds may be ejected by random turbulent motions rather than accelerated by the wind. 
Finally, I show that energy balance cannot be achieved at all for turbulent media above a surface 
density of sa 10^ Mq pc~^. 

Subject headings: galaxies: starburst — ISM: jets and outflows — ISM: structure 



1. INTRODUCTION 

Galaxy outflows are well observed over many masses 
and redshifts [e.g. Martin 1999; Pettini et al. 2001; 
Veilleux, Cecil & Bland-Hawthorn 2005) and form in all 
galaxies in which the star formation rate per unit area, 
E*, exceeds the "Heckman limit" of E^"* « O.IM© yr"! 
kpc~^ (Heckman 2002). They shape the galaxy mass- 
metallicity relation {e.g. Tremonti et al. 2004), affect the 
gas content and number density of dwarf galaxies {e.g. 
Scannapieco et al. 2001; Benson et al. 2003), and enrich 
the intergalactic medium {e.g. Scannapieco et al. 2006). 
Yet, they are extremely difficult to model theoretically, 
as a complete model of galaxy outflows requires a deep 
understanding of both star formation and the evolution 
of the multiphase interstellar medium (ISM) . 

Despite these complexities, I show in this Letter that 
many of the features of galaxy outflows can be under- 
stood from a more general perspective. Regardless of 
how stars impact their environment, a generic result of 
their evolution is the turbulent stirring of the ISM, and 
the resulting turbulent support is at least as important as 
thermal or magnetic pressure (Elmegreen & Scalo 2004). 
Unlike other sources of support however, turbulence re- 
quires both continuous energy input, as it decays away 
to thermal energy, and continuous radiative cooling, as 
this thermal energy must be removed from the system. 
Here I show that this equilibrium is stable only when 
E* < O.IM0 yr~^ kpc~^, while galaxies with larger E^, 
values will develop outflows with many of the key fea- 
tures observed in optical and X-ray studies. 

The structure of this work is as follows. In §2, I show 
how E^"^'* arises from a turbulent heating instability that 
occurs due to the decreasing efficiency of radiative cool- 
ing above 10^'^K. In §3, I explore the implications of this 
instability for the properties of both starbursting galaxies 
and the hot and cold phases observed in galaxy outflows 
(where "cold" means < lO^K in this context). Finally, in 



§4, I describe an extreme case in which turbulent heat- 
ing cannot be sustained for a single dynamical time, and 
relate it to the maximum observed stellar density. 

2. CRITICAL SURFACE DENSITY 

In this section, I show that many of the features of out- 
flowing galaxies can be understood as arising from three 
underlying properties: hydrostatic equilibrium, turbu- 
lent support, and self-regulation. Each of these concepts 
can be formalized quantitatively. 

Hydrostatic equilibrium relates pressure and density 
as ^ w « ^TT/OgGEo, where p and pg are the pressure 
and gas density in the midplane of a disk (or center of 
a spherical distribution), H is the gas scale height, Eq 
is the total matter surface density, and G is the gravita- 
tional constant. This can be rewritten as a function of 
the gas surface density. Eg as 



p w 47rGE„Eo. 



(1) 



Turbulent support can be formalized by requiring that 
turbulence provides a fixed and substantial fraction / of 
the total pressure: 



fp ^ Pgalo, 



(2) 



where gid is the average ID turbulent velocity disper- 
sion. In a galaxy, auj can often be supersonic (Mac Low 
& Klessen 2004), such that / w 1, but it cannot be ex- 
tremely subsonic, such that / ^ 1 {e.g. Elmegreen & 
Scalo 2004; Dalcanton & Stilp 2010). 

Finally, in a disk galaxy one can formalize the concept 
of self-regulation with the Toomre parameter, Q, which 
is above 1 if the pressure is sufficient to stabilize pertur- 
bations on scales too small to be stabilized by rotation. 
In a pure gas disk, Qgas = cri£)K(/^/^7rGEg)~-'^, where k 
as is the epicyclic frequency (Toomre 1964), while in a 
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disk of gas and stars 
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where R — a^j ctid is the ratio of stellar and gas velocity 
dispersions, = Eq — Sg is the stellar surface density, 
and the denominator is evaluated at its maximum as a 
function of wavenumber fc, with q = kaiu /k. (Jog & 
Solomon 1984; Rafikov 2001). In nearby disk galaxies 
Sh./So > 1 and i? > 1, or S^/Sq < 1 and i? > 10 (Leroy 
et al. 2008; Romeo & Wiegert 2011). Over this full range, 
the denominator above varies from w 0.7 — 1, becoming 
very close to 1 if i? w 1 or ^ Sg, as is often the case in 
starbursts. Thus I approximate Qgas+stars as QgasSg/So 
and set Qgas+stars 1, to give: 



(4) 



The motivation for this requirement is both theoret- 
ical, because Qgas+stars ^5 1 leads to stirring of turbu- 
lence by star formation and spiral density waves {e.g. 
Binney & Tremaine 2008) and observational, because 
Qgas+stars ~ 1 — 2 is measured over a wide range of 
disk galaxies (Leroy et al. 2008). From eqs. (JlJ, ([2]), 
and (HI) , this distribution will have a vertical scale length 
H « (Tid/(4/^/^k). For spherical systems, an analogous 
requirement can be derived directly from hydrostatic 
equilibrium as Pgcr^jj/f ~ ^■nGYioTig implies gid ~ 
/i/^TrGSo^ff and H w CTi£,%/(4/i/2)^ where H and 
tg = y^A/irGpo are the radius and the free-fall time. 

Together these three properties lead to a further con- 
straint. Turbulence is continuously cascading to smaller 
scales and thermalizing on a eddy turnover time, Lt/ a^Y>, 

where ctsd = S^/^ctid and Lt and is the length scale of 
the largest eddies. This means that turbulence must be 
continuously driven, and furthermore that thermal en- 
ergy must be continuously removed from the system. In 
a galaxy, this removal occurs through radiative cooling, 
which can be described by a cooling function A that de- 
pends on both gas temperature T and metallicity Z. This 
gives A{T,Z){pg/fimpf « pga^j^Lr'^ > pgoj^H-^, 
where p,mp is the average particle mass and the inequal- 
ity arises because the largest turbulent eddies must be 
smaller than the scale height. Rearranging this in terms 
of the Eg gives: 



?D, (5) 

a constraint that is completely independent of how tur- 
bulence is driven. 

Furthermore, the cooling function cannot grow without 
bound, but rather it is divided into two regimes. Below 
10^-^K, A increases strongly with temperature, reaching 
a maximum of K^^^ « 10~^^ erg cm^ s~^ for material 
with Z K, Zq. In this range, gas cooling is relatively 
stable, and a density perturbation in pressure equilib- 
rium will cool at a similar rate as its surroundings. On 
the other hand, above 10^-^K, a density perturbation 
will cool much quicker than the surrounding gas. This 
regime is unstable and cannot be maintained indefinitely. 
Rather, turbulence will amplify density contrasts, lead- 
ing to a distribution of cold clumps imbedded within a 



hot diffuse medium (Mathews & Bregman 1978; Fall & 
Rees 1985). 

In this case, within the hot regions, cooling becomes 
increasingly less efficient because both density drops and 
A(T, Z) decreases with increasing temperature. The re- 
sult is that the temperature grows to several times the 
post-shock temperature for a v = a^D shock, increasing 
the sound speed to well above the escape velocity and 
leading to a global outflow. This turbulent runaway has 
been simulated in Scannapieco, Gray, & Pan (2012; here- 
after SGP12) and it is fundamentally different from the 
picture of a hot superbubble that bursts out from within 
a cold medium {e.g. Mac Low & Ferrara 1999; Silich & 
Tenorio-Tagle 2001). Instead, hot gas is continuously 
generated over many dynamical times and (although the 
numerical diffusion in SGP12 was too high to demon- 
strate this conclusively) the hot and cold gas phases are 
likely to exist cospatially at all times, as different aspects 
of the same multiphase, unstable turbulent medium. The 
10^'^K dividing line at which turbulent runaway drives an 
outflow corresponds to (Ti£) w 35 km/s (SGP12). From 
eq. this immediately gives a critical surface density 
above which outflows will be generated of 



Y^outflow 
^0 



300kio /-i/^Mopc"^ Disk 
300tg\o/"i/2MQpc-2 Spherical, 



(6) 



where kiq = lOMyrs x k and tff.io = iff/10 Myrs. 

For disk galaxies, Kennicutt (1998) found that = 

2.5 ± 0.7 X lO-^Mgyr-ikpc"^ (Eg/lMg pc-2)^-^=^°-^^ 

and Ei — O.OITE^^, where fl is the circular ve- 
locity. For a flat rotation curve one can use these 

2/5 

fits to obtain kiq = 1.3 (Eg/IOOM© pc~^) , and a 
fit to the data of Leroy et al. (2008), Table 7, gives 
1 (Eo/100Mqpc-2)-i/4 



Sg/Eo 



11±5 



near E™***°™, such 



that Eg"*fl°" 



= (200 ± 50)/-5/7Mqpc" 



3. IMPLICATIONS 
3.1. Implications for Starhursting Galaxies 

The generation of galaxy outflows by unstable tur- 
bulent support immediately leads to several important 
implications. Assuming that thermal pressure, turbu- 
lence, and magnetic fields are in rough equipartition, 
/ « 1/3, we have Eg > 55 ± 35M0 pc"^ and E* > 
0.08 ± 0.06Moyr-ikpc"^ which is in excellent agree- 
ment with E"'* Ki O.lAf0 yr"-'^ kpc^^. This limit corre- 
sponds to a minimum pressure in an outflowing star- 
burst of nT w 7 X lO^Kcm"'^, and unlike models that 
compare the energy input to the gravitational binding 
energy (e.g. Dekel & Silk 1986), it does not depend on 
galaxy mass, which explains the presence of outflows in 
massive Lyman-break galaxies (Shapley et al. 2003), 
luminous and ultraluminous infrared galaxies (Martin 
2005), and blue-cloud star- forming galaxies (Weiner et 
al. 2009). Only two previous papers have attempted to 
explain this limit: Shu, Mo, & Mao (2005) who associate 
it with a velocity of 160 km/s below which they claim 
winds are undetectable (see however Martin 2005), and 
Murray, Menard, & Thompson (2011), who associate it 
with the value at which giant molecular clouds disrupt at 
velocities above the galaxy escape velocity. This, unlike 



3 



our model, predicts a strong scaling of S"'* with galaxy 
mass. 

Secondly, our model predicts that the turbulent veloc- 
ity scales roughly with Eq, and should be at least 35 km/s 
in all starbursts. Recently, Genzel et al. (2011) compiled 
a large sample of the higher-redshift data and found that 
all galaxies with > S"'* do in fact have ciD.Ha > 35 
km/s. 

Thirdly, in our model, the ISM of all outflowing galax- 
ies contains gas with a wide range of temperatures in 
rough pressure equilibrium, including gas above a few 
million degrees, the post-shock temperature of gas col- 
liding at a few times the mean critical turbulent velocity 
(SGP12). Again, these features are well supported ob- 
servationally. For example, recent observations of M82 
uncover a medium in which w 10^ K dense clumps 
are surrounded by > lO^'K gas. In both components, 
nT w 2 X 10"^ K'cm-3 (Strickland & Heckman 2007). 
Note that this vahic, which is similar to that seen in other 
starbursting galaxies (Heckman et al. 1990), is a few 
times greater than the minimum pressure for turbulent 
runaway, suggesting that if SGP12 had included super- 
nova material (e.g. Hill et al. 2012) in addition to turbu- 
lence, even hotter gas may have added to the hot medium 
initially opened up by turbulent runaway (Strickland & 
Heckman 2009). 

3.2. Implications for Galaxy Outflows 

For most galaxies, this hot gas will move many scale 
heights per Myr, while a typical starburst continues for 
many Myrs {e.g. Greggio et al. 1998; Forster Schreiber 
et al. 2003). Thus the hot gas distribution should be well 
approximated by an equilibrium configuration. Such a 
solution was first derived for a constant energy and mass 
input rate within a spherical distribution (or equivalently 
for a conical distribution with a fixed opening angle) by 
Chevalier and Clcgg (1985), providing model that is often 
applied for simple estimates. 

This can be significantly improved as described in 
Silich et al. (2011), by accounting for the profile of the 
turbulent gas distribution that adds energy and gas to 
this wind. In this case the equations of mass, momen- 
tum, and energy conservation are: 



^-^iPhVhr )=qn 
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(8) 



(9) 



where p;, , v;,, P;,, and 7 are the density, radial veloc- 
ity, and ratio of specific heats of the hot gas. The mass 
and energy input per unit volume are given by truncated 
exponential profiles: 



9m = 



_ / gm,oexp(-r/rc) 



0, 



if r < rmax 



(10) 



and (7e,o = 9m,oCQ f^/ (7 — 1), where co,/i > 500 km/s is the 
sound speed of the hot gas at r = 0, and Vc, r^ax are the 
scale radius and truncation radius, respectively. 



These equations can be easily solved, yielding the mod- 
els shown in Figure 1. In the top row of this figure, I plot 
the hot gas velocity, Vh, and sound speed Ch as a function 
of radius for models with rmax = 2rc through 8rc- For 
small Tmax values, the wind properties change dramati- 
cally at rmax, while for large rmax values, the radial trend 



is more gradual. In the 



= 2rc model, for example. 



the sound speed within the driving region is roughly con- 
stant, but drops off quickly outside of the sonic point at 
r = rmax- In models with larger truncation radii, Ch drops 
more significantly within rmax and the transition at rmax 
becomes less pronounced. Finally, when rmax ^ "ire the 
sonic point stays fixed at « Arc, rather than moving out 

to rmax- 

The center row of this figure shows the normalized den- 
sity profile of the hot gas Ph{r) = ph{r)/ph{0)- Unlike the 
sound speed and velocity, ph (r) shows a dramatic change 
at rmax for all values of rmax/^'c- In fact, the model with 
the largest truncation radiiis shows a sharp change in 
phir) at rmax, but no feature at all at the sonic point. 
Note that increasing the sound speed at r = would 
have no effect on ph{r), and would raisc^ both Cs{r) and 
Vh{r) by a fixed factor without changing their radial de- 
pendence. 

Finally, the bottom panel of Figure 1 estimates the 
cold gas ejection from these galaxies. While the cold 
gas found in galaxy outflows is thought to be accelerated 
by the hot wind {e.g. Martin 2005), there is a serious 
problem with this picture. Klein et al. (1994) studied 
the evolution of nonradiative cold clouds propagating at 
a velocity v through a hot medium, showing that if v is 
much greater than the cloud's internal sound speed, it 
will be shredded on a "cloud crushing" timescale 



tc. 



X 



cloud 



(11) 



where x = Pcioud/p/i is the density contrast between the 
media, and Rdoud is the cloud radius. Subsequent studies 
showed that magnetic fields {e.g. Mac Low et al. 1994) 
and radiative cooling (Fragile et al. 2005) can only delay 
this disruption by 1-2 cloud crushing times. On the other 
hand, to accelerate the cloud to the hot wind velocity, the 
impinging mass must be comparable to its total mass, 
which takes a time 



accel 



(47r/3)i?gi„„dPcioud 



t:R: 



loudPhV 



(12) 



Because they arc in rough pressure equilibrium (Strick- 
land & Heckman 2007), the density contrast between the 
< lO^'K clouds and the hot gas is x ~ Ih/T'coid > 1000, 
such that t accel exceeds 30 cloud crushing times. Given 
the results discussed above, it would appear that cold 
clouds will never survive to reach velocities comparable 
to that of the hot wind. 

There is, however, one likely caveat. If cloud cooling 
is efficient and the relative velocity exceeds the hot gas 
sound speed, a bow shock develops in front of the cloud. 
This protects the cloud from ablation, both by reduc- 
ing heating within the cloud and shear around its sides. 
When combined with cooling, these mitigating effects al- 
low the cloud to remain intact for many cloud crushing 
times. In fact, the acceleration of clouds without strong 
disruption in this case has been simulated in both Cooper 
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Fig. 1. — Equilibrium wind solutions with truncated exponential mass and energy input profiles. From left to right the columns 
correspond to models with truncation radii, rmax = 2rc, 4rc, and 8rc, respectively, labeled by their approximate escape velocities Vcac ~ 
'''ID \/2''max//»'c, ~ 85 km/s rmax/'^c where I assume (jid ~ 35 km/s and / = 1/3. Top row: Radial velocity (solid) and sound speed 
(dashed) of the hot gas as a function of radius. The shaded regions show the gas within the sonic point, which lies at the truncation radius 
if ''max $ 47*0 and at lies Arc for larger rmax values. Center row: Normalized hot gas density profiles. Bottom row: Normalized density 
of the cold gas reservoir (solid), and the cold gas expected to be ejected by turbulent motions (dashed). 



et al. (2009) and Kwak et al. (2011), although these runs 
were terminated before the cloud velocity approached Vh ■ 
Nevertheless, we can infer that the acceleration of cold 
clouds is likely to be efficient if they enter the hot gas 
outside of the sonic point, and inefficient if they enter 
the flow where Vh < Ch- 

In Figure 1, we see that only galaxies with escape ve- 
locities > 200 km/s have a portion of their cold gas reser- 
voir outside of the subsonic region. For these, one would 
observe many cold gas clouds with velocities comparable 
to Vh outside of the sonic point, but significant shred- 
ding of cold clouds inside of the sonic point. In fact, 
the archetypal high- mass starburst, M82, displays ex- 
actly these features. Within 300 pc of the disk. X-ray 
observations indicate significant mass loading of the hot 
wind (Suchkov et al. 1994) and the cold gas motions are 
minimal (McKeith et al. 1995). Outside 300 pc, the cold 
gas motions quickly approach Vh, and strong broad line 
SII emission is observed (Westmoquette et al. 2009), indi- 
cating a large population of cold clouds interacting with 
the hot outflow. 

Similarly, the fastest cold gas velocity as measured by 
Nal absorption approaches Vh for galaxies with circular 
velocities > 150 km/s (Martin 2005), while the fastest 
cold gas velocity is ^100 km/s in lower mass galax- 
ies (Martin 2005). These correspond to models with 
?'max < 4rc, for which there is no cold reservoir beyond 
the sonic point. However, these potential wells are small 
enough that the stochastic turbulent motions themselves 
will accelerate a fraction of the gas to above the escape 
velocity. The lower panel of Figure 1 shows the proflle of 



cold gas expected to be ejected by turbulence, calculated 
simply as jOeject(?') = pdr) 0.5 erfc[wcsc 
I crudely approximate v, 



{r)/y/2al^], where 



(0)^1 -J^Amax. For 

these galaxies, it is easy to understand the often perplex- 
ing observation that the hot gas is likely to be moving 
at > 500 km/s while the ejected cold clouds are moving 
< 100 km/s {e.g. Martin 2005). These are young clouds 
ejected by ISM turbulence, and they have not yet expe- 
rienced signiflcant acceleration by the wind. If they had, 
they would have already been shredded. 

4. EXTREME SURFACE DENSITY LIMIT 

Finally, there is an extreme regime, in which Eg is sim- 
ply too high to be stable to the Toomre criterion with 
signiflcant turbulent support, while at the same time bal- 
ancing energy input with cooling, even if A = A™'''^. In 
this case eqs. (U) and ([5]) cannot be simultaneously satis- 
fled for a single dynamical time, and there is no way for 
the gas to achieve even an unstable equilibrium. For disk 
and spherical configurations this occurs when the surface 
density is above 



vmax 
^0 



3x 105 (Eg/Eo)i/2MQpc- 
3 X 10^ (Eg/Eo)i/2M0pc- 



Disk 

Spherical, 
(13) 

where I have taken / = 1/3. 

These values closely correspond to the maximum ob- 
served stellar surface density of Emax ~ lO^M^Qpc"^, 
which is consistent across globular clusters, nuclear star 
clusters, super star clusters, ultra-compact dwarfs, and 
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the centers of elliptical galaxies (Hopkins et al. 2010). 
Above Sg"^''' S^-s can only be driven out within a dynam- 
ical time or fall freely without achieving equilibrium. In 
fact, Larson (2009) showed that if gas above this density 
were deposited onto a central supermassive black hole, 
its mass would match the observed relation to the stellar 
velocity dispersion. Perhaps, like galaxy outflows, this 
is a consequence of the physics of turbulent support at 
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